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A  wide  range  of  experimental  data  on  the  oxidation  of  ZrB2  and 
HfB2  as  a  function  of  temperature  (800°-2500°C)  is  interpreted 
using  a  mechanistic  model  that  relaxes  two  significant  assump¬ 
tions  made  in  prior  work.  First,  inclusion  of  the  effect  of  volume 
change  associated  with  monoclinic  to  tetragonal  phase  change  of 
the  Me02  phases  is  found  to  rationalize  the  observations  by 
several  investigators  of  abrupt  changes  in  weight  gain,  recession, 
and  oxygen  consumed,  as  the  temperature  is  raised  through 
the  transformation  temperatures  for  Zr02  and  Hf02.  Second, 
the  inclusion  of  oxygen  permeability  in  Zr02  is  found  to  ratio¬ 
nalize  the  enhancement  in  oxidation  behavior  at  very  high  tem¬ 
peratures  (>1800°C)  of  ZrB2,  while  the  effect  of  oxygen 
permeability  in  Hf02  is  negligible.  Based  on  these  consider¬ 
ations,  the  significant  advantage  of  HfB2  over  ZrB2  is  credited 
to  the  higher  transformation  temperature  and  lower  oxygen 
permeability  of  Hf02  compared  with  Zr02. 


Symbol 
A  W/ A  (kg/m2) 
W0JA  (kg/m2) 
JT(m) 

L  (m) 

h  (m) 


hex t  (m) 

T  (K) 
t(s) 

Mt  (kg/mol) 

Pi  (kg/m3) 

/ 

4) 


Nomenclature 

net  change  in  weight  per  unit  area 
weight  of  02  consumed  per  unit  area 
recession  of  MeB2 

thickness  of  zirconia  region  in  the 
scale 

thickness  of  zirconia  scale  over  which 

B203(7)  is  present 

external  thickness  of  B203(7) 

temperature 

time 

molecular  weight  of  species  i 
density  of  species  i 

effective  fraction  of  pores  in  Me02 
that  is  permeable  to  gas 
pore  volume  fraction 
tortuosity  factor 


N.  S.  Jacobson — contributing  editor 


Manuscript  No.  25581.  Received  December  2,  2008;  approved  February  9,  2009. 

This  work  was  supported  in  part  by  USAF  Contract  #  FA8650-04-D-5233. 
tAuthor  to  whom  correspondence  should  be  addressed,  e-mail:  triplicane.parthasarathy@ 
wpafb.af.mil 


rP  (m) 
r\  i  (Pa  s) 

J,  (mol  (m2  •  s)  '  ) 

]J  (mol  (m  s  -  atm)-1) 


a 

Pi  (atm) 

T  (mol  (m2  •  s)-1) 

Avap,,  (mol  (m2  •  s)-1) 

Dt  (m2/s) 

£>1,(2, ...  ,0  (m2/s) 

-Do,,Me02  (m2/s) 

A  (m2/s) 

Ct  (mol/m3) 

9 

F  (A  •  s/mol) 

g0  (Q-1  •  m-1) 

g,-  (Q_1  •  m_1) 

Superscripts 

a 

zb 

s 


pore  radius 
viscosity  of  species  i 
flux  of  species  i 

permeability  coefficient  of  oxygen  in 
liquid  boria 

effective  evaporation  coefficient 
partial  pressure  of  species  i 
rate  of  addition  of  boria  to  the 
external  boria  scale 
rate  of  evaporation  of  species  i  at  the 
external  surface 
diffusivity  of  species  i 
diffusivity  of  species  1  in  a  multicom¬ 
ponent  gas  mixture 
diffusivity  of  oxygen  in  Me02 
Knudsen  diffusivity 
concentration  of  species  i 
h/L 

Faraday  constant 
oxygen  ionic  conductivity 
hole  conductivity 

ambient 

zirconia-boria  interface 
MeB2-Me02  interface 
B203(/)-B203(g)  interface 


I.  Introduction 

Diborides  of  Zr  and  Hf  are  now  well  recognized  as  the  most 
promising  refractory  materials  that  are  resistant  to  envi¬ 
ronmental  degradation  in  air  at  very  high  temperatures,  espe¬ 
cially  with  SiC  additions.1^  However,  little  or  no  work  has  been 
reported  on  quantitative  models  of  oxidation  for  this  class  of 
materials.  A  few  recent  works  focused  on  thermodynamic  stud¬ 
ies,1,5  which  have  identified  the  key  condensed  and  gaseous  ox¬ 
idation  products.  Modeling  the  kinetics  of  oxidation  for  this 
class  of  materials  has  been  the  focus  of  our  research  program. 
Experimental  data  on  the  kinetics  of  oxidation  of  SiC-contain- 
ing  diborides  are  only  beginning  to  appear  in  the  literature, 
while  a  large  body  of  data  exist  on  the  pure  diborides.6  11  Ex¬ 
istence  of  this  large  data  base  combined  with  the  simpler  scale 
morphology,  make  modeling  the  oxidation  of  the  pure  dibo¬ 
rides,  the  logical  first  step  toward  understanding  this  class  of 
materials. 
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From  an  initial  study,  we  reported  on  a  model  for  the  oxi¬ 
dation  of  these  materials  in  the  intermediate  temperature  re¬ 
gime,  viz.  ~  1000°-1800°C.12  We  later  extended  this  to  other 
temperatures.13  By  comparing  the  models  with  the  experimental 
data,  we  showed  that  the  data  can  be  interpreted  by  a  model  that 
assumes  the  formation  of  a  porous  zirconia,  whose  pores  tend  to 
be  filled  with  liquid  boria (/),  with  the  boria (/)  evaporating  as 
B203(g)  at  higher  temperatures. 


ZrB2  +  (5/2)02  -<■  ZrO 2{s)  +  B203(/) 

T  >- 

-  873  K 

B2o3(/)  -  b2o3(s) 

T  >- 

-  1273 K 

ZrB2  +  (5/2)02  -  ZrO 2(j)  +  B203(g) 

T  >- 

-  2073  K 

Below  1273  K,  the  model  predicted  an  external  boria  to  form 
and  it  showed  that  inclusion  of  its  flow  was  important  to  ratio¬ 
nalize  data  at  temperatures  above  ~  1073  K.  Above  1273  K,  the 
boria  evaporation  was  predicted  to  be  significant  resulting  in  a 
porous  zirconia  that  is  partially  filled  with  liquid  boria.  Above 
~2073  K,  the  boria  evaporation  was  predicted  to  be  faster  than 
the  rate  of  formation  of  boria,  resulting  in  the  scale  being  es¬ 
sentially  porous  zirconia.  This  model  helped  rationalize  a  large 
portion  of  the  data  reported  in  the  literature.  However,  a  few 
important  effects  were  not  captured.  In  particular,  it  is  found 
that  these  works  could  not  rationalize  some  of  the  abrupt  tran¬ 
sitions  observed  in  the  high  temperature  data,  as  reported  in  the 
literature.  These  transitions  occur  around  1400  K  for  ZrB2,  and 
around  2000  K  for  HfB2.  Second  the  recession  rates  reported  by 
Clougherty  in  air  at  temperatures  above  ~2273  K  were  found  to 
be  significantly  larger  than  was  predicted  by  the  model  for  ZrB2. 
Third,  the  model  could  not  simultaneously  predict  both  the 
weight  gain  and  oxygen  consumed  reported  for  low  tempera¬ 
tures  (1273-1473  K)  for  ZrB2  by  Tripp  and  Graham. 

In  this  work,  we  retain  all  elements  of  the  prior  model,  but 
relax  two  of  the  assumptions.  First,  we  show  that  by  including 
the  effect  of  the  phase  transformation  of  the  Me02  oxide  scale 
from  monoclinic  to  tetragonal  with  accompanying  volume 
change,  the  abrupt  transitions  observed  in  reported  experimen¬ 
tal  data  can  be  predicted/rationalized  accurately.  Inclusion  of 
this  mechanism  is  also  found  to  be  important  in  obtaining  good 
predictions  of  both  the  weight  gain  and  oxygen  consumed  si¬ 
multaneously  for  ZrB2.  Including  the  permeation  of  oxygen 
through  the  Me02  oxide  is  shown  to  rationalize  the  acceleration 
in  oxidation  at  the  very  high  temperatures  for  ZrB2,  while  this 
oxygen  flux  contribution  is  shown  to  be  negligible  for  the  case 
of  HfB2. 

The  experimental  data  on  these  materials  come  from  a  variety 
of  sources  that  used  different  methods  to  characterize  the  kinet¬ 
ics,  viz.,  weight  gain,  interface  recession  and/or,  oxygen  con¬ 
sumed.  The  materials  have  also  been  studied  under  a  variety  of 
different  environmental  conditions  and  temperature  regimes. 
The  model  is  able  to  rationalize  all  of  these  data. 


II.  The  Model 

The  key  elements  of  the  model  described  in  detail  in  our  earlier 
work  are  retained.12  As  a  brief  overview,  a  schematic  of  the 
general  model  is  shown  in  Fig.  1,  which  illustrates  the  assumed 
morphology  of  the  oxidation  products,  representing  steady-state 
conditions.  The  assumptions  are  based  on  experimental  obser¬ 
vations  as  noted  in  our  prior  work.12  The  oxidation  product, 
viz.,  the  scale,  consists  of  Me02  and  liquid  boria,  with  B203(g) 
as  the  dominant  gaseous  species.  At  low  temperatures  where  the 
vapor  pressure  of  B203  is  very  low,  the  pores  between  the  col¬ 
umnar  Me02  product  are  filled  with  liquid  boria,  which  is  re¬ 
tained  by  the  low-surface  tension  and  good  wettability  of  Me02 
by  boria.  At  low  temperatures  the  scale  supports  an  external 
boria  film.  This  external  boria  that  forms  at  low  temperatures 
can  flow  under  gravity  or  due  to  shear  forces  from  external  fluid 
flow.  Oxygen  diffuses  through  the  external  boria  layer  and  then 
through  the  tortuous  pathway  between  zirconia  grains  to  reach 
the  substrate. 


B203(l)  b2o3(!) 


ZrB2 

ZrCV/] 

Fig.  1.  Schematic  sketch  of  the  oxidation  products  and  morphology 
assumed  in  the  model.  The  progression  of  the  cartoons  represents  steady 
state  for  increasing  temperature. 


As  the  temperature  is  raised,  volatilization  as  B203(g)  results 
in  loss  of  the  external  boria  film,  so  that  the  porous  channels  of 
the  Me02  scale  are  then  partially  filled  with  liquid  boria.  Gas¬ 
eous  oxygen  must  then  diffuse  through  the  outer  part  of  the 
zirconia  pore  channel,  and  then  permeate  through  the  boria  in 
the  base  of  pore  channel  to  reach  the  substrate.  The  model  as¬ 
sumes  that  a  quasi-steady  state  is  reached  where  the  rate  of  boria 
formation  is  balanced  by  the  rate  of  boria  evaporation,  which 
determines  the  depth  to  which  the  zirconia  pores  are  filled  with 
boria.  At  yet  higher  temperatures,  where  boria  evaporates  as 
fast  as  it  forms  so  that  liquid  boria  is  absent  in  the  pores,  ox¬ 
idation  is  limited  by  gaseous  diffusion  of  oxygen  in  the  pores, 
which  has  a  weak  temperature  dependence.  In  addition,  the 
evaporation  of  Me02  becomes  significant  resulting  in  a  reduc¬ 
tion  in  the  oxide  scale  thickness.  This  enhances  the  oxidation 
rate,  and  is  therefore  taken  in  to  account. 

Thus  the  rate-limiting  step  for  oxidation  at  intermediate  tem¬ 
perature  is  the  permeation  of  oxygen  through  liquid  boria  in 
series  with  gaseous  diffusion  of  oxygen  through  porous  open 
channels.  Oxygen  permeation  through  solid  Me02  is  negligible 
at  low  and  intermediate  temperatures,  due  to  the  low  electronic 
conductivity  of  Zr02  and  Hf02.  Diffusion  in  the  pores  is  as¬ 
sumed  to  follow  Fick’s  law,  and  kinetic  laws  are  derived  for 
planar  symmetry. 

In  this  work,  we  retain  derivations  from  our  prior  model  for 
the  low-temperature  regime,  but  relax  two  of  the  assumptions, 
one  in  the  intermediate  temperature  regime  (~  1273-2073  K) 
and  another  in  the  high-temperature  regime  (above  ~2073  K). 
In  the  intermediate  temperature  regime,  we  rationalize  the  effect 
of  the  phase  change  of  Me02  from  monoclinic  to  tetragonal  to 
effect  a  change  in  the  effective  porosity  of  the  scale  (Fig.  2).  The 
transformation  results  in  volumetric  shrinkage  with  increasing 
temperature.  We  suggest  that  the  monoclinic  Me02  nucleates 
and  grows  during  the  transient  heat  up  of  the  sample;  as  the 
sample  crosses  the  transformation  temperature  the  Me02  in  the 
scale  transforms  to  tetragonal,  shrinking  in  the  process  and  ex¬ 
panding  the  gap  between  the  Me02  grains.  Thus  we  introduce 
an  increase  in  porosity  of  the  scale  at  the  transformation  tem¬ 
perature,  and  we  assume  that  this  higher  porosity  is  retained 
during  further  oxidation.  Second,  at  very  high  temperatures 
when  the  scale  is  free  of  boria  (typically  above  ~2073  K),  the 
permeation  of  oxygen  ions  through  the  Me02  is  considered  and 
included  in  the  total  flux  of  oxygen  through  the  oxidation  scale. 


(1)  Model  Equations 

The  present  work  extends  that  detailed  in  Parthasarathy 
et  al .  13  The  reader  is  referred  to  those  works  for  detailed  der¬ 

ivations.  Here  the  derivations  are  limited  to  what  is  new;  but  for 
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Fig.  2.  A  schematic  sketch  showing  the  model  assumptions.  During 
transient  heat-up  the  oxide  scale  consists  of  m- Zr02.  As  the  temperature 
rises  above  the  transformation  temperature,  the  m- Zr02  converts  to 
t- Zr02  with  an  associated  shrinkage,  which  opens  up  the  pathway 
for  oxygen  permeation,  enhancing  oxidation.  The  model  assumes  that 
this  opening  remains  constant  as  oxidation  continues  to  thicken  the 
scale. 


continuity,  equations  derived  in  our  prior  work12,13  and  used  in 
this  work  are  briefly  summarized.  All  the  model  variables  are 
listed  and  defined  in  the  nomenclature. 

The  following  shows  equations  for  the  three  different  regimes. 
In  all  regimes,  the  diffusivity  of  gases  in  a  multicomponent  gas 
mixture  and  the  effect  of  Knudsen  diffusion  were  modeled  as 
detailed  in  Parthasarathy  et  al.12  At  the  outer  surface,  either 
boria  or  Me02  evaporation  takes  place  under  boundary-layer- 
limited  diffusion.  This  evaporation  at  the  surface  of  species  i, 
with  molecular  weight  Mh  at  partial  pressure  Pt,  at  temperature 
T,  was  taken  to  follow  the  modified  Langmuir  law  that  includes 
the  effect  of  boundary  layer  as  follows: 


•W  =  u.Pi/(2nMjRT)0'5 


(2) 
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where  Cdopant  is  the  trivalent  dopant  (impurity)  concentra¬ 
tion  in  the  Me02,  aQ  is  the  ionic  conductivity,  and  cj^(Po2)  is 
the  hole  conductivity  of  Zr02  at  partial  pressure  of  oxygen, 
P0l.  The  oxygen  diffusivity  in  hafnia  has  been  reported  to 
be  a  factor  of  about  5  lower  than  that  of  zirconia.15  Thus  the 
ionic  conductivity  of  hafnia  was  taken  to  be  one -fifth  that  of 
zirconia. 

Noting  that  the  ionic  conductivity  is  known  to  be  Po2  inde¬ 
pendent  and  hole  conductivity  has  a  (Po2)1^  dependence,  the 
flux  of  oxygen  through  the  oxide,  Me02,  is  now  given  by  Wag¬ 
ner’s  equation  as,16 


|^02,Me02 


1  -f  RT  f2 
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^^dlnPo2 
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where  Pq7  and  PsQi  are  the  oxygen  partial  pressures  at  the  am¬ 
bient  and  at  the  oxide/boride  interface,  and  f  is  the  fraction  of 
pores  in  the  oxide.  Combining  Eqs.  (3),  (4),  and  (5)  and  inte¬ 
grating,  one  obtains 


\Jo2. 


Me02 
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where  aQ,  cj^(Pq  )  are  given  by  Eqs.  (3)  and  (4).  The  parallel 
flux  of  oxygen  through  the  gas  phase  in  the  pores  is  given  as 


|/o 


2  gas 


=  fD0l 


(7) 


The  net  flux  of  oxygen  to  the  MeB2/Me02  interface,  and 
its  relation  to  flux  of  B203  and  formation  rate  of  Me02,  is  given 
by 


(|7o  2,Me02  |  +  | /o  2gas|)  ^  —  |«Zb203  I  —  ^Me02 


(8) 


where  a  represents  the  effect  of  boundary-layer-limited  diffu¬ 
sion,  and  varies  with  experimental  conditions,  but  is  reported  to 
be  in  the  range  of  10-3-10-2  under  most  laboratory  condi¬ 
tions.14  A  value  of  5  x  10-3  was  used  for  all  model  results  shown 
in  this  work. 

(A)  High-Temperature  Regime:  As  pointed  out  in  Part¬ 
hasarathy  et  al.}2  the  permeabilities  of  oxygen  through  zir¬ 
conia  and  hafnia  at  high-oxygen  activities  are  dominated  by 
electron-hole  conductivity,  which  is  much  lower  than  the 
ionic  conductivity  except  at  very  high  temperatures.  Thus  as¬ 
suming  that  Me02  is  inert  to  oxygen  permeation  in  the  inter¬ 
mediate  temperature  regime,  as  carried  out  in  Parthasarathy 
et  al.}2  was  a  good  assumption.  However,  at  very  high  temper¬ 
atures,  above  ~2000  K,  the  hole  conductivity  becomes  signifi¬ 
cant  and  oxygen  permeation  through  Me02  must  be  taken  into 
account.  In  this  work,  the  model  for  the  high-temperature  re¬ 
gime  presented  in  Parthasarathy  et  al.13  is  modified  to  include 
this  effect. 

The  effect  of  diffusion  of  oxygen  through  the  zirconia  phase 
in  the  scale  was  taken  into  account  as  follows.  As  elaborated  in 
Parthasarathy  et  al.}2  for  Zr02,  the  ionic  and  hole  conductiv¬ 
ities  were  taken  as 


where 
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where  Z>o2,  Db2o3  are  the  gas  phase  diffusivities  of  oxygen  and 
B203  through  the  pores,  and  Reaction  is  the  chemical  equilibrium 
constant  for  the  reaction,  MeB2+5/202->Me02+B203(g). 

Equation  (7)  along  with  Eqs.  (5,  6)  and  (8,  9)  can  be  solved 
numerically  for  Pb2o3  and  ^o2 .  These  values  were  then  used  in 
the  following  differential  equation,  which  was  numerically  inte¬ 
grated  to  obtain  scale  thickness  L ,  recession  X ,  weight-gain  per 
unit  area  AW/ A,  and  weight  of  oxygen  consumed  per  unit  area 
WoJA. 
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Note  that  the  model  calculates  the  contribution  of  ambipolar 
oxygen  diffusion  in  the  Me02  phases  assuming  that  the  oxide 
pegs  retain  a  dense  morphology  of  constant  cross  section.  Thus 
the  calculated  oxygen  flux  contribution  should  be  considered  the 
maximum  possible  flux  contribution. 

(B)  Intermediate  Temperature  Regime:  The  phase  trans¬ 
formation  of  zirconia  is  assumed  to  result  in  a  change  in  the 
effective  volume  fraction  of  pores  in  the  zirconia.  Because  the 
transformation  from  monoclinic  to  tetragonal  results  in  volu¬ 
metric  shrinkage,  the  fraction  of  pores  will  increase  as  the  tem¬ 
perature  reaches  the  transformation  temperature.  Thus  in  the 
model  we  have 


/  =/i  for  T  <  rtrans 
f  =fi  for  T  >  Ttrans 


with /2</i.  The  other  parameters  are  given  by  equations  derived 
in  Parthasarathy  et  al.}2  summarized  below,  for  the  key  param¬ 
eters  of  interest,  the  internal  depth  of  liquid  boria,  h ,  the  zirconia 
scale  thickness,  L,  recession  of  MeB2,  and  weight-gain  per  unit 
area,  A  W/A,  and  weight  of  oxygen  consumed  per  unit  area, 
WqJA 
h  =  qL ; 
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(C)  Low-Temperature  Regime:  The  parameters  of  inter¬ 
est  are  external  boria  thickness,  the  thickness  of  the  zirconia 
scale,  recession,  and  weight  gain.  In  the  absence  of  viscous  flow 
of  the  external  boria,  the  effective  external  layer  thickness  is  de¬ 
termined  by  the  rate  of  production  of  boria  less  the  rate  of 
evaporation  of  boria.  The  boria  included  in  the  underlying  po¬ 
rous  zirconia  layer  also  contributes  to  the  effective  diffusion 
barrier  thickness. 
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When  viscous  flow  is  significant,  the  external  boria  layer 
thickness  is  decreased  to  a  value  determined  by  the  rates  at 
which  it  flows  and  thins  at  the  surface,  as  given  below.  Assuming 
laminar  flow  of  the  boria,  the  “effective”  external  thickness,  hex t, 
is  given  by  Perry  and  Chilton17 


hex  t  — 


3  (Mb2q3  T  b2Q3  /Spec)flB2Q3 
gpl2o,  sin(<l>) 
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where  TB2o3  is  the  number  of  moles  of  boria  formed  and  added 
per  unit  area,  per  unit  time,  on  the  surface. 
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The  Me02  scale  thickness,  recession  of  MeB2,  weight-gain  of 
the  sample,  and  the  weight  of  oxygen  consumed  are  then  given 
by 
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III.  Model  Predictions  and  Validation 

(1)  Parameters 

A  list  of  variables  used  in  the  model  as  shown  in  the  nomencla¬ 
ture  with  a  brief  description.  The  model  uses  data  from  the  lit¬ 
erature  for  all  thermodynamic  quantities,  which  are  available  in 
the  compendium  by  Barin.18  The  diffusivities  of  gases  in  mul¬ 
tigas  mixtures  were  calculated  using  parameters  given  by 
Svehla,19  and  the  methods  outlined  in  Parthasarathy  et  al.12 
The  Knudsen  effect  on  diffusivity  was  included  as  determined  by 
the  pore  radius.  The  temperature-dependent  viscosity  of  boria 
was  obtained  from  the  works  of  Eppler20  and  Li  et  al.21  The 
permeability  of  oxygen  in  liquid  boria  was  obtained  from  the 
works  of  Tokuda  et  al?2  and  Luthra.23 

The  only  unknown  parameters  in  the  model  are  the  pore  ra¬ 
dius,  effective  (tortuosity-included)  volume  fraction  of  pores  be¬ 
fore  and  after  the  phase  transition  of  Me02,  and  the 
concentration  of  dopants  (unintentional  impurities)  in  the 
Me02.  The  pore  radius  was  taken  to  be  0.5  pm,  which  is  based 
on  microstructures  reported  in  the  literature;  further  parametric 
studies  in  Parthasarathy  et  al.12  showed  that  there  is  no  signifi¬ 
cant  effect  unless  the  pores  are  <0.1  pm  in  size.  The  trivalent 
dopant  concentration,  C^opant,  in  Me02  was  taken  to  be  0.1%, 
which  is  reasonable  given  the  purity  levels  of  the  diborides  used 
in  the  experimental  work.  For  the  case  of  Zr02,  at  the  highest 
temperatures  of  the  test,  raising  the  dopant  levels  to  0.5%  gave 
better  fit  to  experimental  data,  implying  impurities  gained  from 
the  furnace  atmosphere  as  a  possible  way  to  rationalize  the  data. 
The  choice  of  dopant  concentration  had  no  effect  on  predictions 
for  HfB2  oxidation,  implying  that  the  effect  of  oxygen  through 
the  Hf02  is  negligible. 

The  effective  volume  fraction,/,  of  pores  is  determined  by  the 
actual  volume  fraction  of  pores,  4>,  and  the  tortuosity  of  the 
pores,  x.  The  tortuosity  is  defined  as24,25 


z  = 
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D 

D^f'1 


D  = 


(iO  ; 


(28) 


Simulations  show  that  the  tortuosity  factor  is  about  25  for 
cj>  =  0.05,  and  ~4  for  <}>  =  0.15.26  Thus  the  effective  volume 
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fraction,  f  can  be  much  smaller  than  the  actual  porosity  at  low 
volume  fractions  of  porosity.  If  the  tortuosity  is  taken  in  to  ac¬ 
count,  values  for  the  effective  volume  fraction,  f  of  0.03-0.04 
(representing  §  ~0.15),  are  very  reasonable,  compared  with  ex¬ 
perimental  images  of  the  porous  scale. 

The  effect  of  phase  transformation  of  Me02  was  estimated 
from  published  data  on  volume  change  associated  with  mono¬ 
clinic  to  tetragonal  phase  transformations.  The  temperature  of 
transformation  was  taken  to  be  1400  K  for  Zr02  and  2000  K  for 
Hf02,  in  close  agreement  with  the  findings  of  the  review  pub¬ 
lished  by  Wang  et  al.f1  who  quote  an  average  value  of  1367  and 
2052  K,  respectively.  The  volume  change  associated  with  the 
transformation  has  been  reported  to  be  ~  3%,  with  the  value  for 
Zr02  being  slightly  larger  than  that  for  Hf02,  in  the  review  by 
Taylor.28  The  effective  pore  fraction,/ that  gave  the  best  fit  was 
0.04  for  Zr02  and  0.03  for  Hf02,  above  the  transformation 
temperature.  Below  the  transformation  temperature,  the  values 
assumed  were  0.0015  and  0.003  for  Zr02  and  Hf02,  respectively. 
The  smaller  pore  fraction  is  attributed  to  the  lack  of  phase 
transformation  for  the  monoclinic  phase,  and  the  finite  value  of 
pore  fraction  is  attributed  to  the  low  surface  tension  of  boria  (at 
600°C,  it  is  0.07  N/m,29  same  as  water  at  25°C)  that  is  likely  to 
result  in  complete  wetting  of  the  zirconia/hafnia  by  liquid  boria. 


(2)  Comparison  of  Model  Predictions  with  Experimental 
Data 

There  are  extensive  data  in  the  literature  on  the  low-tempera¬ 
ture,  intermediate  temperature,  and  high-temperature  regimes  of 
oxidation  of  ZrB2  and  HfB2,  although  there  are  a  limited  num¬ 


ber  of  investigators.  Berkowitz-Mattuck10  has  conducted  exten¬ 
sive  tests  and  reported  on  the  temperature  dependence  of 
parabolic  rate  constants  for  recession  and  oxygen  consumed  us¬ 
ing  gas  conductivity  measurements  under  37.5  torr  (0.0049 
MPa)  partial  pressure  of  oxygen  in  helium  carrier  gas.  Weight 
gain  and  oxygen  consumed  as  a  function  of  time  in  the  low- 
temperature  regime,  under  250  torr  (0.033  MPa)  of  pure  oxygen, 
have  been  reported  by  Tripp  and  Graham.11  Interface  recession 
and  oxygen  consumed  at  higher  oxygen  partial  pressures  (air)  in 
the  intermediate  and  high-temperature  regimes  have  been  re¬ 
ported  by  Kaufman  et  al?°  For  the  high-temperature  oxidation, 
the  recession  rates  in  1  h  in  air  were  measured  by  Kaufman  and 
Clougherty.6  8  Kaufman30  has  also  reported  on  the  time  depen¬ 
dence  of  oxygen  consumed  at  high  temperatures  in  0.005  MPa 
02  partial  pressure  for  ZrB2.  We  examine  the  validity  of  the 
model  by  comparing  with  these  data  sets. 

Figure  3  shows  the  data  for  ZrB2  and  HfB2,  along  with  model 
predictions  for  recession  rate  constant  and  the  parabolic  rate 
constant  for  oxygen  consumed.  The  jumps  in  the  data  are  well 
represented  by  the  model,  thus  attributing  these  transitions  to 
the  volumetric  shrinkage  upon  phase  change  for  Zr02  and 
Hf02,  At  higher  temperatures,  where  there  is  no  liquid  boria 
in  the  scale,  the  model  includes  the  diffusion  of  oxygen  ions 
through  the  oxide  in  addition  to  gaseous  diffusion  through  the 
pores.  As  seen  from  Eqs.  (2)  and  (3)  the  permeation  flux  is  de¬ 
pendent  on  the  concentration  of  the  dopant  (here  impurities)  in 
the  oxide.  For  ZrB2,  the  fit  with  the  experimental  data  at  the 
highest  temperatures  is  good  with  0.5%  impurity  in  the  oxide, 
indicating  that  the  experimental  samples  had  some  impurity  (not 
surprising)  or  that  the  samples  picked  up  impurity  from  the  en- 


1 /Temperature,  K 


1 /Temperature,  K 


Fig.  3.  The  effect  of  phase  transition  of  the  Me02,  and  the  effect  of  oxygen  permeation  through  Me02  at  very  high  temperatures  are  shown  by 
comparing  model  predictions  to  experimental  data.  The  plots  compare  the  model  predictions  of  parabolic  constants  for  recession  and  oxygen  consumed 
as  a  function  of  temperature  for  ZrB2  and  HfB2,  with  experimental  data  reported  in  the  literature.  The  model  predicts  the  transitions  based  on  the  phase 
transition  of  Zr02  at  1400  K  and  that  of  Hf02  at  2000  K,  which  is  very  close  to  phase  transition  temperatures  reported  in  the  literature.  At  the  highest 
temperatures,  using  0.5%  dopant  in  Zr02  resulted  in  a  better  fit,  while  permeation  of  oxygen  through  Hf02  was  found  to  be  negligible  even  at  1  %  dopant 
concentration. 
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Fig.  4.  Model  predictions  for  time-dependent  behavior  are  compared  with  reported  data.  For  the  low- temperature  regime,  the  weight  gain  and  oxygen 
consumed  as  a  function  of  time,  reported  on  ZrB2  in  0.033  MPa  oxygen,  by  Tripp  and  Graham  are  shown.  For  the  intermediate  regime,  the  oxygen 
consumed  for  a  sample  tested  in  0.0026  MPa  oxygen  helium  mixture  at  1711  K  is  compared  with  the  data  of  Berkowitz-Mattuck.  For  the  high- 
temperature  regime,  the  data  for  HfB2  in  0.005  MPa  Q2/He  mixture  of  Kaufman  is  compared  with  the  model  prediction. 


vironment  during  high-temperature  oxidation  treatments.  The 
value  chosen  for  impurity  level  had  no  effect  on  the  predictions 
for  HfB2,  and  the  fit  with  the  experimental  data  is  good,  imply¬ 
ing  that  HfB2  is  insensitive  to  impurities. 

Figure  4  shows  the  model  predictions  for  time-dependent  ox¬ 
idation  behavior  in  the  low-temperature  regime,  compared  with 
those  measured  by  Tripp  and  Graham11  on  ZrB2  samples  under 
0.033  MPa  oxygen  partial  pressure  of  pure  oxygen.  Both  the 
weight  gain  and  the  oxygen  consumed  are  shown  to  be  in  good 
agreement  indicating  that  the  model  interpretations  are  reason¬ 
able.  The  data  for  intermediate  temperature  are  due  to  Berko¬ 
witz-Mattuck,10  and  the  high- temperature  data  on  HfB2  is  from 
Kaufman  et  al.30  The  model  is  shown  to  be  very  reasonable  un¬ 
der  all  temperature  regimes. 

Figure  5  shows  the  dependence  on  oxygen  partial  pressure 
predicted  by  the  model  compared  with  the  data  of  Berkowitz- 
Mattuck  for  HfB2.  The  oxygen  partial  pressure  dependence 
shows  a  transition  from  nearly  independent  at  low  oxygen  par¬ 
tial  pressures  to  nearly  linear  at  higher  partial  pressures  of  ox¬ 


ygen.  The  mechanism  for  this  transition  arises  from  the 
transition  in  oxygen  transport  being  dominated  by  permeation 
through  liquid  boria  at  lower  partial  pressures,  to  diffusion  of 
gaseous  oxygen  through  pores  in  the  zirconia  at  higher  partial 
pressures.  Thus  the  transition  in  the  dependence  with  respect  to 
Po2  or  temperature,  is  related  to  the  fraction  of  the  oxygen 
diffusion  path  along  liquid  boria  versus  gaseous  pore.  Higher 
oxygen  partial  pressure  and  lower  temperatures  favor  retention 
of  liquid  boria  and  thus  a  linear  dependence  on  oxygen  partial 
pressure,  as  was  shown  in  Parthasarathy  et  al.12 


IV.  Discussion 

This  work  examined  the  effects  of  the  shrinkage  of  the  oxide 
scale  due  to  phase  transformation  and  of  the  higher  temperature 
ambipolar  diffusivity  of  oxygen  in  the  oxide  scale  on  the  kinetics 
of  oxidation  of  ZrB2  and  HfB2. 
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Fig.  5.  The  effect  of  oxygen  partial  pressure  on  the  measured  parabolic  rate  constant  of  oxygen  consumed  is  compared  with  the  model  predictions  for 
HfB2,  at  two  different  temperatures. 


The  model  introduced  the  effect  of  phase  transformation  us¬ 
ing  a  simple  adjustment  to  the  pore  fraction  in  the  oxide  scale 
formed,  but  in  proportion  and  magnitude  to  experimentally 
observed  volume  change  and  the  actual  temperature  of  the 
transformation.  The  agreement  with  experimental  data  are 
found  to  be  very  good.  This  suggests  that  the  effect  of  phase 
transformation  of  the  Me02  from  monoclinic  to  tetragonal 
plays  an  important  role  in  the  oxidation  kinetics  of  MeB2.  The 
much  better  resistance  to  oxidation  for  HfB2  in  the  temperature 
of  experimental  measurements  (typically  <1600°C,  the  trans¬ 
formation  temperature  of  Hf02)  can  be  attributed  entirely  to 
the  higher  transformation  temperature  of  2000  K  for  Hf02 
versus  1400  K  for  Zr02.  If  this  is  true,  then  additives  that  in¬ 
crease  the  transformation  temperature  or  suppress  the  trans¬ 
formation  temperature  must  be  beneficial  to  the  oxidation 
resistance  of  ZrB2  and  HfB2. 

The  model  introduces  the  effect  of  ambipolar  diffusion  of 
oxygen  in  Me02  at  higher  temperatures  using  the  semiempirical 
expressions  derived  in  Parthasarathy  et  al.12  This  requires  that  a 
small  percentage  of  extraneous  dopant/impurity  of  a  different 
valency  be  present  in  the  oxide  and  that  its  amount  in  the  scale 
be  known.  Assuming  a  value  of  0.1%  of  trivalent  dopants  in  the 
oxide  scale  is  able  to  rationalize  the  high-temperature  data  well. 
For  ZrB2,  a  value  of  0.5%  rationalizes  data  at  the  highest  tem¬ 
peratures  tested,  implying  the  possibility  of  impurities  picked  up 
during  oxidation.  At  2273  K,  the  relative  contribution  of  ambi¬ 
polar  diffusion  of  oxygen  to  the  total  oxygen  permeation  is 
about  30%  at  0.1%  impurity  level  and  rises  to  about  60%  at 
0.5%  impurity  level.  It  is  noteworthy  that  HfB2  is  insensitive  to 
impurity  level  in  its  oxide;  the  relative  contribution  of  ambipolar 
diffusion  to  total  oxygen  ingress  is  <  1  %  at  2273  K.  Once  again, 
the  better  oxidation  resistance  of  HfB2  over  ZrB2  suggested  by 
the  model  is  credited  to  the  lower  diffusivity  of  oxygen  in  hafnia 
compared  with  zirconia,  as  reported  by  Smith  et  al}5  In  addi¬ 
tion,  the  lower  vapor  pressures  of  hafnia,  minimizes  evaporative 
loss  in  the  extreme  high-temperature  regime. 

In  summary,  a  model  that  includes  phase  transformation  of 
the  oxide  scale  and  the  ambipolar  diffusion  of  oxygen  through 
the  oxide  at  high  temperatures  is  found  to  rationalize  the  ex¬ 
perimental  data  reasonably  well.  Alloying  additions  that  sup¬ 
press  or  delay  the  phase  transformation  would  be  beneficial. 


V.  Summary 

A  wide  range  of  experimental  data  on  the  oxidation  of  ZrB2  and 
HfB2  as  a  function  of  temperature  from  800°  to  2500°C  are  in¬ 
terpreted  using  a  mechanistic  model  that  relaxes  two  assump¬ 
tions  made  in  prior  work.  First,  inclusion  of  the  effect  of  volume 
change  associated  with  monoclinic  to  tetragonal  phase  change 
of  the  Me02  is  found  to  rationalize  the  observations  by  several 
investigators  of  abrupt  changes  in  weight  gain,  recession  and 
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oxygen  consumed,  as  temperature  is  raised  through  the  trans¬ 
formation  temperature  in  both  ZrB2  and  HfB2.  Second,  the  in¬ 
clusion  of  oxygen  permeability  in  Zr02,  neglected  in  prior 
treatments,  is  found  to  help  rationalize  the  behavior  at  very 
high  temperatures  (>  1800°C)  of  ZrB2,  where  an  enhancement 
in  oxidation  has  been  observed.  At  2273  K,  the  relative  contri¬ 
bution  of  ambipolar  diffusion  in  Zr02  to  oxygen  ingress  is  about 
30%  at  1%  impurity  level  in  the  oxide.  The  effect  of  oxygen 
permeability  in  Hf02  is  negligible  (<  1  %)  and  thus  the  oxidation 
of  HfB2  is  insensitive  to  trivalent  impurities.  Based  on  this  anal¬ 
ysis,  the  significant  advantage  of  HfB2  over  ZrB2  is  credited  to 
the  higher  transformation  temperature  and  lower  oxygen  per¬ 
meability  of  Hf02  compared  with  Zr02. 
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